Random unitary dynamics of quantum networks 
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We investigate the asymptotic dynamics of quantum networks under repeated applications of random unitary 
operations. It is shown that in the asymptotic limit of large numbers of iterations this dynamics is generally 
governed by a typically low dimensional attractor space. This space is determined completely by the unitary 
operations involved and it is independent of the probabilities with which these unitary operations are applied. 
Based on this general feature analytical results are presented for the asymptotic dynamics of arbitrarily large 
cyclic qubit networks whose nodes are coupled by randomly applied controlled-NOT operations. 
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In recent years specific features of classical networks have 
been investigated intensively fit] because they exhibit generic 
organizing principles shared by rather different systems, such 
as living cells or the internet. In view of these developments 
and in view of recent significant progress in the control and 
manipulation of quantum systems |2(| it is natural to extend 
these investigations into the quantum domain and to explore 
characteristic phenomena of quantum networks. 

In a natural generalization of its classical analog the nodes 
of a typical quantum network are formed by spatially local- 
ized distinguishable quantum systems and couplings between 
different nodes originating from interactions or communica- 
tion are described in general by a completely positive quan- 
tum operation 01 ■ Iterative applications of this quantum op- 
eration give rise to a dynamical evolution of this quantum net- 
work. So, in general contrary to typical interacting many- 
body quantum systems f% the couplings between different 
parts of a quantum network cannot be described by infinites- 
imal generators, such as Hamiltonians or Lindblad operators. 
Quantum networks can model rather different physical quan- 
tum systems, such as interacting gases or in the context of 
quantum information processing a quantum-communication- 
based internet. Contrary to their classical counterparts they 
are not only able to manipulate classical correlations but also 
to distribute entanglement which is a characteristic quantum 
resource. As strong entanglement cannot be shared freely be- 
tween several quantum systems iU the distribution of entan- 
glement over a quantum network follows rules which are fun- 
damentally different from those governing the distribution of 
classical correlations. 

In this letter we explore characteristic features of the in- 
tricate interplay between entanglement and decoherence in 
quantum networks. Whereas unitary evolution typically gen- 
erates entanglement, i.e. genuine quantum correlations, de- 
coherence tends to destroy them. Thus the competition be- 
tween these two counteracting tendencies is expected to pro- 
duce interesting dynamical features. In order to explore some 
of them let us consider the rather simple case of quantum 
networks whose nodes are formed by elementary two-level 
quantum systems (qubits) coupled by randomly chosen uni- 
tary controlled-NOT operations. In large quantum networks 
with large numbers of coupled nodes this rather elementary 
random unitary dynamics gives already rise to complicated 



dynamical features under repeated iterations. Even for mod- 
erately large numbers of qubits (between 10 and 30) a detailed 
description becomes computationally intractable. Therefore, 
alternative ways have to be developed. 

In the following it is shown that despite these complications 
the asymptotic dynamics of iterated random unitary quantum 
operations exhibits very regular patterns. It turns out that their 
asymptotic dynamics is governed by a typically low dimen- 
sional attractor space which is determined completely by the 
unitary transformations involved and which is independent of 
the probability distributions with which these unitary transfor- 
mations are selected. This basic property of random unitary 
quantum operations allows us to investigate systematically the 
asymptotic dynamics even of large quantum networks. As an 
example we present first results exploring characteristic fea- 
tures of cyclic qubit networks whose dynamics is determined 
by random unitary controlled-NOT operations. 

In order to put the problem into perspective consider an ar- 
bitrary quantum system whose dynamics is described by the 
iterated application of a random unitary operation 

Thus, a single application of this quantum operation consists 
of selecting a unitary linear operator Ui = {JJi \ I E 1} 
randomly according to the normalized probability distribu- 
tion {0 < pi < 1, i G 1} and of applying it onto an ini- 
tially prepared quantum state p G B{Ti). {B{Ti) denotes 
the Hilbert space of linear operators over the d-dimensional 
Hilbert space Ti. describing the quantum system under con- 
sideration.) Correspondingly, after n iterations of this ran- 
dom unitary operation the quantum system is in the state 
= $($($(... ((l>(p))))). 

The dynamics of Eq.([T]i can model various physical situa- 
tions. For example, if the quantum system consists of a large 
number of distinguishable few-level systems it can model col- 
liding distinguishable low-dimensional systems whose colli- 
sions are well separated in time. In this case each Ui e U 
describes the unitary evolution of a completed collision which 
occurs with probability pi. Alternatively the iterated dynam- 
ics of Eq.([T]i can also describe a quantum walk |6i] whose ideal 
unitary dynamics is perturbed by random imperfections. In 
the context of quantum information processing this dynam- 
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ics may model the probabilistic exchange of quantum infor- 
mation between different nodes within a quantum internet in 
which different nodes entangle each other by unitary opera- 
tions selected randomly by its users. 

Random unitary operations have a number of characteristic 
properties. Obviously they belong to the class of unital quan- 
tum operations iHQl which leave the maximally mixed quan- 
tum state p = invariant. Defining the Hilbert-Schmidt 
scalar product [A, B) — Ty{A^ B) in the Hilbert space of lin- 
ear operators B{T-l) the adjoint of $ of Eq.([T]i is given by 

^Hp) = Y^P^UlpU,. (2) 

Thus, in general $ is not normal, i.e. [<1>^,$] 7^ 0, so that 
it cannot be diagonalized. Nevertheless, its properties can 
still be analyzed systematically with the help of Jordan nor- 
mal forms 1811 . We show in the following that despite the 
resulting complications the asymptotic dynamics of iterated 
random unitary operations is governed by surprisingly 
regular patterns which considerably simplify their description 
especially in the large-n limit. In special cases this asymptotic 
dynamics can even be determined analytically. 

For a discussion of this asymptotic dynamics we use the 
fact that <I> of Eq.([T]i cannot decrease the von Neuman entropy 
S, i.e. 

Si^p)) > Y.p..s{u,pUl)^Sip). (3) 

The resulting monotony of S ($(")(/5)) and its boundedness 
for the finite dimensional quantum systems at hand imply that 
in the hmit of large numbers of iterations (1>'^"' (p) leads to a 
constant von Neuman entropy. Thus, the relations 

tj^p^Uj = U,,pnUl (4) 

have to be fulfilled for all asymptotic quantum states /5„ = 
^'^"^ (p) with n sufficiently large and for all unitary operations 
with i, i' E I. From Eq.® we note that the set of all its solu- 
tions forms a linear space, the attractor space A C B and that 
$ of Eq.([TJ acts on this attractor space unitarily. As unitary 
transformations are normal the restriction of $ onto the attrac- 
tor space A can be diagonalized with the help of a complete 
set of orthonormal eigenvectors Xx e A. These eigenvectors 
fulfill the eigenvalue equations and orthonormality constraints 

U^XxUj = XX X, Tr (xt^v) = Sxx' (5) 

for all i E I. From Eq.Q it is apparent that the linear op- 
erators of the attractor space A also form a C* -algebra. This 
property implies a number of useful relations, such as 

I A I = 1, xl^Xx*, A" A"" 7^ 1 — > Tr (x^'X;?) = 0. 

Thus, the attractor space A is defined by all linear oper- 
ators Xx fulfilling Eqs.(|5]l with the additional requirement 
I A 1= 1. Furthermore, all asymptotically accessible quan- 
tum states pn = <i>^"-'(p) are contained in this attractor space. 



i.e. lim„_,oo Pn S A. Let us now introduce the (continuous) 
projection operator 

Vi.) ^ Tr(^l {■))Xx (6) 

|A| = 1 

which projects the linear space B onto the attractor space A. 
In view of the relation 1/A* = A for | A| = 1 it commutes with 
the random unitary operations of Eq.([T]i, i.e. [P, $] = 0, so 
that we can conclude with the help of Eqs.Q and (|6]l 

lim ip)^v( lim (p)) = lim (p^ut) = 
lim [//;'p„„«J7^;)t ^ ^ ^n^r (Xl {p) ) Xx (7) 

|A| = 1 

with Iq G I arbitrary and with the projected initially prepared 
quantum state jjout = 'P{p)- Eq.© is a main result of this 
letter It shows explicitly how the asymptotic dynamics of $ 
of Eq.([T]i is determined by the structure of the attractor space 
A. Of course, Eq.Q could also have been derived in a direct 
but more cumbersome wa y w ith the help of a Jordan normal 
form decomposition of $ lUOll . 

According to Eq.(|7]i the determination of the asymptotic dy- 
namics of can be divided into two main steps. In the first 
step one determines the set of eigenvalues A of unit modu- 
lus and the associated orthonormal basis {Xx] of the linear 
space A with the help of Eqs.Q. Typically, this is a difficult 
task which may be facilitated by symmetries. In the second 
step one evaluates the unitary asymptotic action of onto 
the projected quantum state pout according to Eq.Q. Typi- 
cally, the attractor space is expected to be of sufficiently low 
dimensions so that Eq.(|7|i implies significant simplifications. 

Eq.(|7J hints at some remarkable general features of the 
asymptotic dynamics of First of all, the space of attrac- 
tors A is determined completely by the set of unitary operators 
U. In particular, this implies that the asymptotic dynamics is 
independent of the classical probability distribution character- 
izing $. Nevertheless, in general this probability distribution 
still influences the rate of convergence. Furthermore, if one 
of the unitary operations from the set U is the unit-operation 
(apart from a global phase) the only possible eigenvalue is 
given by A = 1. The resulting asymptotic dynamics is thus 
stationary. Another general consequence can be derived for 
unitary operations U which generate a ray-representation of a 
finite group. Because the unit-element of this group is a prod- 
uct of the generators contained in U, for each eigenvalue A 
there is a natural number nx with A"^ = 1. Thus, the result- 
ing asymptotic dynamics is periodic. However, Eq.(|7J also 
applies to the most general non-stationary and non-periodic 
cases of asymptotic dynamics. 

Let us now apply the general results of Eq.(|7]i to the de- 
scription of the asymptotic dynamics of a particular class 
of quantum networks consisting of qubits which are cou- 
pled by randomly selected controlled-NOT operations. A 
controlled-NOT operation between qubits i and j is defined by 
Ci,j\a)i®\b) j = I a) 1 6© a) j. It is known to entangle qubits 
effectively because it produces pure maximally entangled Bell 
states from separable states of the form (|0)i + |l)i)/ ^/2®\b)j, 
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for example. Here, the pure states \a) with a G {0, 1} denote 
orthonormal basis states of the computational basis of a qubit 
and a © 6 with b G {0, 1} denotes addition modulo 2. Let us 
further assume that our network has a one-dimensional cyclic 
topology so that nodes i and i + 1 with i = 1 , . . . , iV are cou- 
pled by controlled-NOT operations Ci,i+i and that in view of 
the cyclic topology qubits + 1 and 1 are identical. 

In order to determine the asymptotic limit of the corre- 
sponding random unitary operation one has to solve the eigen- 
value problem of Eqs.(|5]l. The possible eigenvalues can be 
determined easily by noting that controlled-NOT operations 
have the property Cf^_^_i = 1. Therefore, the only possible 
eigenvalues of Eqs.© are given by A = ±1. The determina- 
tion of all eigenvectors with eigenvalue A = 1 is facilitated by 
the observation that the pure quantum states |0) and |$) = 

E'"r'k>/V2^^(with \z) ^ M|jiv-i)---b2)|ji>, 
z ~ J2i=i 2*^^ji, and ji G {0, 1}) are invariant under all 
controlled-NOT operations of the cyclic network. Using the 
additional fact that the unit operator is a solution of Eq.© 
with eigenvalue A = 1 we find the following orthonormal 
eigenvectors for A = 1 

Xi = |0)(0| , X2 = |0)($|, ^3 = |$)(0|, (8) 
X4 = |<i>)(<i>| , X5 = (l-|0)(0|-|$)($|)/V2^-2. 

By a somewhat lengthy calculation one can prove from 
Eqs.© by induction that there are no additional eigenvectors. 
Thus, for any number of qubits N the eigenspace associated 
with A = 1 is five dimensional and is given by Eq.®. By 
induction one can also demonstrate that for N > 2 solutions 
of Eqs.© with eigenvalue A = —1 do not exist. It is only 
in the special case of = 2 that a non trivial normalized 
eigenvector exists. Explicitly it is given by ifioll 

X, ^- |0)|1)(1|(0| + |0)|1)(1|(1| + |1)|0)(0|(1|- 

|1)|0)(1|(1|-|1>|1)(0|(1| + |1)|1>(1|(0|. (9) 

Having determined an orthonormal basis of the attractor 
space A, the general form of the asymptotic dynamics of the 
quantum network can be determined with the help of Eq.O. 
Thus, for any number of qubits with N > 2 the projected 
quantum state pout is given by 

Pout P ^ +'y^^P> 2N_2 ' ^ ^ 

This state is stationary because it is invariant under all 
controlled-NOT operations under consideration. Here, the 
projection operator P2 = |0)(0| + |$)($| projects onto the 
two-dimensional subspace 7^2 C H spanned by the pure 
states |0) and |$) which are invariant under all controlled- 
NOT operations under consideration. The probability p = 
Tr{P2/5} measures the overlap of the initially prepared quan- 
tum state p with this invariant subspace. Eq.(fTOli implies 
that any initially prepared quantum state p which is contained 
completely in subspace B{'H2) is not affected by the randomly 
applied controlled-NOT operations. Thus, B{'H2) forms a 
decoherence-free subspace lfTll[T2ll . 
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FIG. 1: Trace-norm distance D between n-th iter- 
ate p„ and the projected quantum state pout for dif- 
ferent pure six-qubit states p — \tp}{^p\: \^) = 
lOOOOOl) (full), IllllU) (dashdot), |101010) (dot). 
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FIG. 2; Trace-norm distance D between n-th iterate pn and the pro- 
jected qua ntum sta te pout for different pure two-qubit states p = 

WW/VMWl- W = 0.3e'^/^|10) +e'^/^|ll) (full), |10) + 
111) (dashdot). 

Let us finally discuss the convergence towards the asymp- 
totic dynamics. A convenient measure for the deviation be- 
tween the quantum state after n iterations /5„ and the corre- 
sponding asymptotic quantum state IJJ^poutU"^ is the Hilbert- 
Schmidt norm or trace distance D between both states with 

- {pn - U^PoutUr^ rPn - KPoutUZ^). 

In Fig[T]numerical results are presented depicting the trace- 
distance D between the states p„ and pout ~ U^^^Pout 

^lo for 

a cyclic quantum network consisting of iV = 6 qubits and 
for various initial quantum states p with a uniform probability 
distribution. It is apparent that convergence is achieved in a 
strictly monotonic way. Nevertheless, the rate of convergence 
depends on the initially prepared quantum state. 

In Fig|2] analogous numerical results are depicted for a 
cyclic quantum network consisting of two qubits and for the 
trace-distance D between the states /5„ and pout- Contrary to 
the previously considered case, now the asymptotic dynamics 
is not always stationary, i.e. typically pout 7^ U{^,PoutUi^J , be- 
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cause according to Eq.(|9|l there is also a non trivial eigenspace 
with eigenvalue A = — 1. Depending on whether the initially 
prepared quantum state p overlaps with this eigenspace or not 
the asymptotic dynamics of /}„ is non-stationary or stationary. 

In conclusion, we have presented a general method which 
allows us to determine the asymptotic dynamics of quantum 
networks under the influence of iterated random unitary op- 
erations. It is based on the determination of the associated 
asymptotic attractor space. This attractor space is independent 
of the probability distribution involved and is typically low 
dimensional thus simplifying the asymptotic dynamical de- 
scription considerably. The results presented for cyclic qubit 



networks demonstrate explicitly how the properties of the at- 
tractor space determine the resulting global entanglement of 
the quantum network. It is expected that the exploration of 
the structure of attractor spaces associated with the dynamics 
of more general networks will shed light also onto other open 
problems of network dynamics in the quantum domain, such 
as the connection between network topology and the resulting 
asymptotic dynamics. 
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